extern doublereal dznrm2_(integer *n, doublecomplex *x, integer *incx);
